Anomalous diffusion with log-periodic modulation in a selected time interval 
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On certain self-similar substrates the time behavior of a random walk is modulated by logarithmic 
periodic oscillations on all time scales. We show that if disorder is introduced in a way that self- 
similarity holds only in average, the modulating oscillations are washed out but subdiffusion remains 
as in the perfect self-similar case. Also, if disorder distribution is appropriately chosen the oscillations 
are localized in a selected time interval. Both the overall random walk exponent and the period of 
the oscillations are analytically obtained and confirmed by Monte Carlo simulations. 
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The basics of anomalous diffusion on fractal substrata 
were established some decades ago [TH5]- However, it has 
been recently realized that, on some self-similar objects, 
the time behavior of a random walk (RW) is modulated 
by logarithmic-periodic oscillations (see, for example, 6- 
I14|). This kind of modulation is very ubiquitous. Exam- 
ples have been observed in biassed diffusion on random 
systems earthquake dynamics [111 [20], escape 

probabilities in chaotic maps [21], processes on random 
quenched and fractal media [2^H^ , diffusion-limited ag- 
gregates [5S] , growth models [27] , and stock markets [55}- 
I5T] . In spite of this diversity, log-periodicity is considered 
evidence of a discrete scale invariance [32j . which results 
from some inherent self-similarity |55] . 

Interestingly, the above described complex behavior 
can be captured by simple models, as is the case of a sin- 
gle particle moving by jumps between nearest-neighbor 
(NN) sites of a one-dimensional lattice. Indeed, the parti- 
cle mean square displacement exhibits a subdiffusive be- 
havior modulated by log-periodic oscillations, provided 
that the hopping rates are distributed in a self-similar 
way |34j . In this work we introduce disorder, by ran- 
domly shuffling those hopping rates, and investigate the 
effect of this modification on the RW dynamics. We give 
theoretical and numerical evidences that, while subdif- 
fusion survives with the same RW exponent as in the 
ordered substrate, the oscillatory modulation is washed 
out. We also show how the oscillations can be localized 
in selected time windows by a proper choice of disorder 
distribution. We believe that our model incorporates the 
essential mechanisms for anomalous diffusion modulated 
by localized oscillations. In this sense, it is a minimal 
model, which can serve for understanding more realistic 
systems. 

Model - We consider a RW on a one-dimensional lat- 
tice, with transitions between NN sites only. The proba- 
bility to hop from site j to site j + 1 (and from site j + 1 to 
site j) per unit time is denoted by k^j\ For readability, 
we schematically represent these rules by a set of barriers 
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of height /i'^^ — c/k^^^ {j e Z), where c is a constant (see 

Fig.g. 

Ordered substrate - We briefly review the main 
properties of the model introduced in Ref. [33] . For more 
details we refer to that paper. The substrate depends on 
two parameters, L (odd integer greater than one) and A 
(real positive), and is built in stages. The result of every 
stage is called a generation (see Fig|l] for L = 3). The 
zeroth-generation corresponds to the situation in which 
all the hopping rates are identical (fc^J) = go, Vj G Z). 
In the first generation, the hopping rate k'^^^ is set to 
?! (< ?o) for every j = pL — (L + l)/2, with p inte- 
ger. All the other hopping rates remain as in the gen- 
eration zero. This iteration continues indefinitely and, 
in general, the generation n is obtained from the gen- 
eration n — 1 by setting fc'-'' = (< gn-i), for every 
j — pL" — (L" + l)/2, with p integer. From now on we 
will refer to the full self-similar substrate obtained after 
an infinite number of iterations as a completely ordered 
substrate. In this limit, the dynamics become full self- 
similar if = 1 + a (n = 0,1,2,...), where 
L)^"^ is the diffusion constant of the nth generation. This 
condition is equivalent to 

^ = ^ + {l + Xy-^XL\ forz = 1,2,3,..., (1) 

which determine the hopping rates for the completely 
ordered substrate. 

It was also demonstrated that, if A'^x{t) is the RW 
mean-square displacement on the completely ordered 
substrate, for a time such that AU^ < (A^a;(t))^/^ < 
j^l^n+i constant of the order of one), every- 

thing will happen as on the nth-generation substrate 
and A'^x{t) = 2D'-^h, for t in that interval. Moreover, 
because of full self-similarity, the sequence of length- 
dependent diffusion coefficients leads to 

A^x{t) = Kt^-'fit), (2) 

where is a constant, — [2 + log(l + A)/logL)]^^ is 
the RW exponent, and f{t) is a log-periodic function of 
period t = L^/'' (i. e., /(tr) = f{t)). This behavior is 
shown qualitatively in Fig. [2] (dots). 
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FIG. 1. Construction of a substrate with a self-similar distri- 
bution of hopping rates (L = 3). Top: The zeroth generation 
(n = 0) with all the hopping rates equal to go- Center: The 
first generation (n = 1), where the hopping rates qi appear. 
The corresponding barriers are separated by a distance L. 
Bottom: The second generation (n = 2) shows the emergence 
of the hopping rate 52 separated by a distance . In the limit 
n — >■ CX3 a full self-similar distribution of hopping rates is ob- 
tained, which corresponds to a completely ordered substrate. 
Numbers below sites indicate position on the lattice. 
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FIG. 2. (Color online) Mean-square displacement as a func- 
tion of time. The full straight lines have slopes of 1, and 
represent normal diffusion, i. e., A'^x = 2D'">t. The dashed 



line represents the global behavior A 



A log-periodic 



modulation is observed on all time scales for a completely or- 
dered substrate (blue dots) and at short times for a partially 
ordered one (red thick full line). 
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FIG. 3. (Color online) Top: Schematic of a completely or- 
dered substrate (L = 3). Numbers bellow barriers indicate 
NN pair label. Bottom: A partial disorder realization of the 
substrate in the top (p = and m = 1). The locations of 
go and qi (black full lines) remain fixed. The others hopping 
rates (red-dashed lines) are distributed at random on the rest 
of the lattice. 



Disorder - For the sake of clarity, when qi is the 
hopping rate between sites j and j + I in the per- 
fectly ordered substrate, we will label pair of sites with 
the number I. Note that the set of NN pairs labeled 
with I can be roughly associated with the length scale 
and that qi appears on the lattice with a frequency 
/, = (L'+i)/(L-l) (see Fig.[|top). 

A disorder realization will implicate a random permu- 
tation of the set {qi\i > 0}, on the same one-dimensional 
lattice. In this context, we will study the disorder av- 
erage of the RW mean-square displacement. It is useful 
to distinguish two different scenarios. When the entire 
array of hopping rates is randomized on the lattice we 
will refer to the resulting substrate as being completely 
disordered. In contrast, we will say that it is partially or- 
dered if only the hopping rates in a subset of the ordered 
array are randomly permuted. 

In each realization of the completely disordered sub- 
strate, the hopping rate between any pair of NN sites 
is qi with a probability Pi = fi (for i > 0), consis- 
tent with the frequency distribution. It is not hard to 
see that the log-periodic modulation, present in the case 
of the self-similar substrate, should be washed out be- 
cause of disorder. For the self-similar case, log-periodicity 
comes from the strong association between the diffusion 
constant D^"^ and length scales between and L"+^ 
(see Fig. |2|, and, thus, should disappear when the spe- 
cial order in the array of hopping rates is lost. On the 
other hand, since we are always dealing with the same set 
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{qi\i > 0}, the average spreading speed (i. e., the speed 
at which a RW explores the lattice) can not be affected 
by disorder Therefore, we expect that the disorder- 
averaged mean-square displacement A'^x behaves, as a 
function of time, as 

A2^(i) = Et^", (3) 

{E is a constant) with the same RW exponent as in the 
perfectly ordered substrate but without oscillations. 

A partially ordered (for length scales between and 
j^p+m+i^ with m > 1) substrate can be defined as fol- 
lows. The hopping rate for every NN pair with its label 
in {p,p + 1, ...jP + m} is the same as in the completely 
ordered case. That is, for any pair labeled with I in that 
set, the corresponding hopping rate is g/. For the other 
NN pairs, qi (with i ^ {p,p + 1, + m}) is randomly 
chosen from the set {qi,withi ^ {p,p+l, ...,p-t-m}} with 
probability P, ^ h/Z, where Z = h- 
Thus, in every realization of disorder, the hopping rates 
associated with length scales between and 2^^+™ re- 
main as in the completely ordered substrate but the other 
are randomly distributed on the rest of the lattice. 

According to the discussion above we expect that 
disorder-averaged mean-square displacement will con- 
tinue to show an overall subdiffusive behavior, with the 
same RW exponent as in the completely ordered sub- 
strate and that, due to partial order, an oscillatory mod- 
ulation appears when BLP < (A^{t)y/^ < BLp+"'+\ 
(where B is a constant). Outside this interval, the oscil- 
lations should be erased by disorder. Note that the geo- 
metrical analysis leading to the relation among L and 
r [34] can be performed at any scale and, thus, r — L^/" , 
also in this case. 

To summarize, we predict a subdiffusive behavior mod- 
ulated by oscillations in a time window determined by 
disorder distribution. For concreteness, let us assume 
p = 0. In this case, the initial ordered array remains 
without changes up to the length scale L™+^. On larger 
scales, the hopping rates are randomly distributed, (see 
the example with L = 3 and m = 1 shown in Fig. |3]). 
We thus expect the behavior sketched in Fig. [2j That is, 
anomalous diffusion (with 2v < 1), modulated by oscil- 
lations at short times and without oscillations at longer 
times. 

Numerical results - To test the predictions outlined 
above, we perform standard MC simulations on a lattice 
of linear size L" {qo = 1/2 and the time step St = 1). 
Every RW starts at the center of the lattice and n is 
always large enough to prevent the RW from reaching 
the lattice borders. 

The disorder-averaged mean-square displacement as a 
function of time, obtained for structures with L = 3 and 
A = 1.48, is plotted in Fig|4j The filled circles correspond 
to complete disorder and the open circles to partial order 
with p — and m — 1. The straight line has a slope 2v 
(with the theoretical value of v) and was drawn to guide 
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FIG. 4. A'^x against t, for a completely disordered substrate 
(filled circles) and for partially ordered one with p — and 
m = 1 (open circles). In both cases, L = 3, and A = 1.48. 
The straight line has the slope 2i/, which represents the over- 
all subdifussive behavior. In the insets we show the scaled 
disorder-averaged mean-square displacement as a function of 
time for each case (symbols as in the main panel). The curve 
in the upper inset corresponds to a harmonic function with 
period r. Both f = 0.354 and r = 22.3 were analytically 
obtained. For more details, see the text. 



the eyes. The good agreement between analytical and nu- 
merical results is apparent in both cases. Let us remark 
that, for the partially ordered substrate, A'^x{t) yields a 
modulated power-law behavior at short times. At longer 
times, the modulation is erased; however, anomalous dif- 
fusion still holds, with the same RW exponent. In the 
insets, we have plotted \og[A'^x/ {Ct'^'^)] against log(t), 
using the same data. The constant C was appropriately 
chosen to have, in the case of partial order, the oscilla- 
tions centered around zero. In the upper inset, we have 
also drawn the curve A sin(27r log(t)/ log(r) -I- a), with 
the theoretical value of t. It is a first-harmonic approxi- 
mation of a log(T)— periodic function, consistent with the 
numerical findings {A and a are fitted parameters). 

We have observed that if the randomness does not af- 
fect the structure on length scales smaller than a cutoff 
distance Ic, the oscillations are present only for times 
shorter than a crossover time tc (~ iV'^)- To get a mod- 
ulation in another time window we explore diffusion a 
substrate with p = 2, m = 1, L = 3, and A = 9. The cor- 
responding numerical results, plotted in Fig[5] show the 
oscillating amplitude localized at the center of the graph; 
outside this region the time dependence of A'^x is well de- 
scribed by a power law. As in previous figures, the inset 
shows a filtered plot and a first-harmonic approximation 
of the data. The modulation is clearly localized between 
two crossover times tci and tc2, with tc2/tci — L^^^"^''/'^ . 
Note that, also in this case, all our predictions are con- 
sistent with the numerical findings (Fig. [s]- inset). 

Conclusions - When disorder is introduced into a 
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FIG. 5. The disorder-averaged mean-square displacement as 
a function of time for L = 3, A = 9, p = 2 and m — 1. The 
straight line has a slope 2u. Inset: A'^x/{Ct^'^) against t for 
the same data. The curve was drawn as in Fig. |4] Here, we 
use tc2/tci = L^^" ■ The RW exponent and the period have 
the analytical values = 0.245 and r — 90.0. 



self-similar structure, a single particle diffusing on the re- 
sulting substrate exhibits an anomalous behavior modu- 
lated by oscillations in some characteristic time intervals. 
Both the RW exponent and the period of the oscillations 
are determined by the parameters of the original per- 
fectly ordered substrate. On the other hand, the position 
and length of the time window in which the oscillatory 
behavior survives randomization depend on disorder dis- 
tribution. Indeed, the modulation holds for those length 
scales that are not affected by disorder; however, it is 
washed out on other scales because of the incoherent su- 
perposition involved in the averaging process. Finally, 
we would like to emphasize that even though we have 
focused on oscillations located in a continuous interval, 
the generalization to the case of several time windows 
separated by gaps of unmodulated behavior is straight- 
forward. Thus, the oscillatory regime can be modified as 
desired by proper choice of L, A and disorder distribution. 
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